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There is, as yet, no complete theory allowing prediction of the 
position of the transition of a laminar compressible boundary layer 
into a turbulent one. Approximated approaches are employed for prac- 
tical purposes. They are generally based on use of the linear theory 
of stability, which has been fairly well developed (see [1] , for 
example) . Use of the linear stability theory in combination with 
information about the field of disturbances provided encouraging 
results for predicting boundary-layer transition on models in wind- 
tunnel tests [2]. In flight tests, however, and in tests in many 
aerodynamics installations , the initial disturbance spectrum in the 
boundary layer is unknown. In this case it is possible to use the 
(rougher, of course) e n method of predicting the transition point, 
which recommends itself well at subsonic speeds in both wind-tunnel 
and flight experiments, including for three-dimensional boundary 
layers (see [3] , for instance) . In this method the position of the 
transition is established upon achievement of disturbance-amplitude 
ratio A = Q/Qq = e (Q 0 is the disturbance amplitude on the lower 
branch of the neutral stability curve, Q is the present value of 
the amplitude) , which is the amplification coefficient of disturbances 
in the unstable region. Since for M > 1 there are no generalizing 
detailed calculations for predicting the transition point, although 
there is a great need for carrying them out, these calculations were 
made . 


Taken as the basis was a program for calculating the coefficients 
for increase in disturbances, a., in a boundary layer with a heat- 

Is 

transfer pressure gradient [4] . The procedure for calculating the 
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stability characteristics is laid out fairly well in [1/ 5] as well. 


The flow of a 
layer is described 


compressible heat-conducting gas in the boundary 
by system of equations 


d (fy u ) , d(r*pt>) n 

ds + dy . 



.dl 


dl 


V u 7s + V V T» = 


d_ (X_ dl \ , d_ f( . X\ d (u 2 ! 2) ' 

dy\ c pdy) d J ^ c p) d y 


(1) 


Here I is the total enthalpy, c is the specific heat, r w ^ s ) is 
the body surface equation, y and X are viscosity and heat conductivity, 
e is a subscript denoting that the parameters are taken at the exter- 
nal border of the boundary layer, k = 1 for the axisymmetric case, 
and k = 0 for two-dimensional flow. From equations (1) , the last 

is valid at dl /ds = 0. 

e 


The flow of air around an impermeable surface with a given wall 
temperature was examined. In this case the boundary additions are /80 

a(0) = iKO) = 0; 1(0) = /«,; H°°) = b(°°) = 1. i 

Assuming limited similitude [6] (as shown by research specially 
conducted by the author, this assumption is quite possible) , system 
(1) was converted into a system of ordinary differential equations 
and numerically integrated (see [1] for more detail). As a result 
of the integration, the distribution of longitudinal velocity and 
temperature, their derivatives, and the viscosity across the boundary 
layer were determined, as required when solving stability equations. 


To determine the coefficients of increase in disturbances, the 
system of stability equations in the Dunn-Lin approximation [7] was 


used : 


p [i (au — ©) u'acp] = — ^ 

VM" Re i 

p [ ia (au — a>) <p] = Srr, 

i (au — a>) r + p'acp + pa (if + q}') = 0, 
p [i (au — to) 9 + aT'cp] = — a (y — 1) (if + q/) 4- 6'. 


( 3 ) 
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(4) 


it t _ j) 

~P ~~P + T 

with boundary conditions 

f(O)=<p(O) = 0(O)=O; 

/, qp, 0 -*■ 0 at y -*■ 

Here u, p, p, and T are the time-averaged velocity, density, pressure, 
and temperature; /, <j> , r, 0, and ir are the disturbance amplitudes of 
longitudinal and normal velocities, density, temperature, and pres- 
sure; a = a + i a. is the wave number of the disturbance in the 

X* i' 

direction of wave propagation; to is the transformed angular velocity 
of the disturbance; Re = Re* cos x, M = Me* cos X, & = a/cos X, 
to = to/cos X/ f = //cos X/ where Re is the Reynolds number, M is the 
Mach number, x is the slope of the wave relative to the main flow. 

The stroke designates differentiation with respect to the coordinate 
normal to the surface. 

The Reynolds number, determined from the transverse (normal) 
scale 6 , 



where u is the longitudinal velocity; p is the density; v is the 
kinematic viscosity (v = y/p) and s is the coordinate along the 
surface of the body; r u (s) is the body surface equation; e is a sub- 
script denoting that the parameters are taken at the external border 
of the boundary layer; k — 1 for the axisymmetric case; k - 0 for 

two-dimensional flow. With gradient-free flow (i.e., at 6 = 

s 2 2k 

(2s /u ) ( du /ds) = 0, where s - fu v p r ds) for the axisymmetric case 

Q & Q & ^ Lc? 

The method of solving system of equations (3) with boundary con- 
ditions (4) is described in detail in [1] . New variables were intro- 
duced, and system (3) was formed into a system of six ordinary first- /81 
order differential equations. It was integrated numerically using the 
orthogonalization method. The calculations assumed Prandtl number 


and for two-dimensional flow. Re 
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o = 0.72, adiabatic constant y = 1.41, and Sutherland's law for 
change of viscosity with temperature. 


As written, problem (3) -(4) corresponds to a two-dimensional 

problem. In fact, however, it allows obtaining information about 

three-dimensional disturbances [7] . The results of integration per- 

2 

mitted obtaining relation a. = $(Re, X/ F ) > where F = to* cos x/ Re ' 

Re = Re/cos X; = oL cos x (X is the angle of disturbance propaga- 
tion) . When solving the problem it was assumed that oo is real, and 

$ ./a . = 3 /a [2], where 3=3 + i3. is the wave number in the lateral 

^ % v v , r> t 

direction. Then x = tan" (3 /a ). 

/V ' ^ / 


The degree of spatial intensification of disturbance is asso- 
ciated with the disturbance amplitude Q by relation 

Real ( 4 M) , ( 6 ) 

From (6) and (5) it follows that 

Ke , 

(* a-dRe 

ln|(?/<? 0 l = -2 J 

2 ^ d(l'A) 


Re. 


( 7 ) 


where 7=1 
I = 1/3. 


For a plate and a wedge, 1=1, for a cone 


Expression (7) defines the ratio of disturbance amplitudes at 
points with coordinates Re and Re Q , which is the coefficient of 
disturbance intensification in the area under consideration. If point 
Re Q corresponds to the lower branch of the neutral stability curve 
and the initial spectrum of boundary- layer natural oscillations is known, 
expression (7) can be used to predict the transition, considering 
(by analogy with M < 1) , for example, that the transition occurs if 
the amount of pulsation at the frequency under examination is on 
the order of 1% of u • 

oo 

If, however, the spectrum of boundary-layer natural oscillations 
is unknown, to predict the position of the transition we can try to 
use the approach to this problem suggested by Smith [8] and developed 
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by a number of authors for the case of subsonic flow. This approach 
is based on the results of work by Schubauer and Skramstad [9] and 
many other works showing that the boundary layer selectively inten- 
sifies the disturbances so that the transition is produced by dis- 
turbances whose frequencies lie in a limited area. For an established 
local Reynolds number, we distinguish from the entire spectrum of 
disturbances observed in the boundary layer the one most magnified 
in amplitude, whose frequency corresponds to the upper branch of 
the neutral curve. Also lying near the upper branch is the frequency 
of the main wave, observed directly in front of the transition. 

Smith's approach is as follows. For an established frequency, 
the maximum of the disturbance-amplitude ratio Q/Qg corresponds to 
the upper branch of the neutral curve. Smith, Jaffe, and Okamura [10] 
suggested associating the position of a transition point found 
experimentally (on various models in low-turbulence wind tunnels) 
with the value of Q/Qg = corresponding to the transition. A num- 

ber of authors have suggested determining the transition from a point 

9 

where the disturbance is magnified e times. At supersonic speeds 
the e n method obviously is also suitable for rough estimation of the 
transition location. Unfortunately, at present it is virtually im- 
possible to check the validity of the choice of specific value for 
n (n = 9, for example) by comparison with experimental data for M > 1, 
because of the absence of full-scale flight test data (except [10] 
at M j< 2) . The value n - 3 [2] is best suited for usual supersonic /82 
wind tunnels (with a fairly high initial level of disturbances in the 
boundary layer) . It must be emphasized once again, however, that 
a more accurate method of determining the transition location should 
include knowledge of the spectrum of external disturbances and laws 
of boundary layer sensitivity (after which the initial spectrum of 
natural oscillations in the boundary layer becomes known) , the ability 
to calculate the nonlinear stage of disturbance development, use of 
a more complete system of stability equations (with all dissipative 
terms, considering non-parallelism of flow) , more accurate determina- 
tion of the criterion of transition beginning (the amplitude "one- 
percent" method of determining the transition is also not entirely 
accurate), and consideration of some other less-significant factors. 
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But for a rough, approximate prediction of the transition point there 
is some sense in using the method employed in this work. 

As opposed to subsonic speeds (where two-dimensional disturbances 
are the most unstable) , at M > 1 we have to examine a. = a. (Re, F, \) , 

"Is 'Is 

where F = to/Re is a dimensionless frequency. The maximum a- for 

"Is 

M = 1.5-4 are reached in the x ~ 50-85° range. In our work the 
calculations were performed for critical angle y* , defined as the angle 

Re 

i'* a ,d Re 
_ 2 \ — 

at which integral ' 1 from (7) most rapidly reached the required 

Reo 

value of n (in particular, n = 9), i.e., at x = X* the Reynolds num- 

Yl 

ber determined at A = e is minimal; it was conditionally adopted as 
the transition Reynolds number Re^. 

The stability characteristics and the conditional location of the 
transition were calculated for the first stable mode of disturbances 
at different Mach numbers (M = 1.5-4) and different values of the 
velocity (pressure) gradient (8 = 0-0.4) on heat-insulated (temperature 

factor T =1) and cooled (T < 1) surfaces. The value of Re. was 

w w t 

calculated chiefly for flat bodies. It should be noted, however, 
that at 8 > 0 calculations were carried out for several arbitrary 
bodies with minor change of 3 and along the surface of the body at 
certain places on it. For a real body with a known distribution of 
3(s) and M^(s) we have to define a conditional location of the transi- 

tion using integral — 2 ] — - As the corresponding calcu- 

lations showed, the difference between the location of the transition 
on an airfoil profile calculated in this manner (from Q/Qq = e ' , 
which is characteristic for wind tunnels) and that obtained in tests 
in a T-313 wind tunnel [13] did not exceed 10%. 

The static temperature for all Mach numbers was taken as the same 
— T = -50°C (save only for calculations intended for comparison of 
theoretical results and experimental data obtained in a wind tunnel) . 
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Results of Calculations 


Figure 1 shows the effect of pressure gradient B on the position 

of the neutral curve (a. = 0) and lines of eaual value of a . (-a . > 0) 

•z- ^ t 

characteristic for a heat-insulated surface. The curves of neutral 

stability are limited below the frequency at which the coefficient 

g 

of amplification A reaches the value of e (this also refers to other 
graphs with coordinates F and Re). The strong stabilizing effect of 
a negative pressure gradient (B > 0) can be seen on both the position 
of the neutral curve (and correspondingly of the critical Reynolds num- 
ber Re cr — the minimum Reynolds number at which disturbances begin 
to increase at any frequency) and the coefficient of disturbance 
increase. Figure 1 also shows the reduction of the frequency de- 
termining the transition with increase in B* 



/83 


Calculations at T =1 and various M and 8 revealed that the 

w 

critical Reynolds number decreases with increasing Mach number in the 
M = 1. 5-4.0 range. A negative pressure gradient (B > 0) (as can be 
seen in Fig. 1) has a strong stabilizing effect on the stability 
characteristics. While for a gradient-free flow the least transition- 
determining frequency corresponds to M = 1.5, at 8 = 0.2-0. 4 it is at 
M = 2.5-4. 
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Based on Fig. 1 and similar graphs for M = 2-4, relations were 
constructed between the conditional transition Reynolds number Re^_ 

and Mach number at T =1, and they are shown in Figs. 2 (3=0) and 

w 

3 (3 = 0.1 and 0.2). Also shown in Fig. 2 are the relations between 
Mach number and critical Reynolds number ^ e cr and maximum coefficients 
of increase in disturbances. A distinct correlation can be seen 
between these two relations and the relation between Re^_ and M. The 
maximum of ratio Re /(-oi^)^^ corresponds approximately to the maxi- 
mum of the relation between Re^_ and M. This is explained by the fact 

that for this case (T =1, 3=0), the lower Re , the lower the Re n 

w cr 0 

Re^ ' i 

in relation — 2 j ajdRe=9 (Re, corresponds to the lower branch of the 

' Re o 

neutral curve, and Re, to the upper) , and the higher -a. in this rela- 

tion, the faster the required value of the integral accumulates, i.e., 

Re t = Re Q + (Re 1 - Re Q ) , where Re Q ^ Re cr and ( Re i “ Re g) ^ V (-^ ) ma x * 

Thus, the lower Re and the greater (-a.) , the lower Re.. And 

cr ^ t max' t 

with increase in 3 at significant rise in Re^., difference 
(Re^ - RCq) - Reg also rises significantly, and at 3 >_ 0.1 the position 
of the transition is determined chiefly by (Re, - Re„) ^ 1/ (-a* ) , 

which is shown by Fig. 3 (Fig. 3a shows the maximum values of the 
coefficients of increase in disturbances at the frequency determining 
the transition) . 



Figure 4 shows a comparison of the values 

of Re^ obtained in this work (at T = 1 and 
t w 

3=0) and experimental data — the results of 
F-15 (cone) flight tests [12] and experiments 
in a T-325 wind tunnel (flat plate) [13] (no 
full-scale flight tests were conducted at 
M > 2) . Calculation of Re for a cone (I = 1/3) 

^ Yl 

(curve 1) were conducted using the e method 

at n = 9. The relation of Reynolds number 

2 7 

determined from A = e (curve 2) is close 
for the results of tests in common supersonic /84 
wind tunnels (where the initial level of dis- 
turbances in the boundary layer is fairly high) 
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a 



(the data in [13] were obtained at a 
stagnation temperature T ^ = 300 K and 

single Reynolds number Re, = ( u/v ) 

6 -L °'' 

= 25 x 10 1/M) . The agreement of experi- 

mental and calculated values of Re^ is 
quite satisfactory, on the whole. It 
should be noted here that exponent 
n =2.7 describes the experimental data 
well at Re x = (25-35) x 10 6 1/M. For 
Re^ x 10 ** = 10 and 60 [1/M] the differ- 
ence between the transition Reynolds 

numbers obtained in the wind tunnel and 

2 7 

calculated from e * reaches 18% at 
M = 4, but in this case the calculated 
and experimental data correlate better 
at another exponent n (2.2 and 3, re- 
spectively) . Both experiments and cal- 
culations demonstrate that there is a 
local maximum (at M - 2) and minimum 
(at M - 3.5-4) in relation Re fc = Re t (M) 

at T , = 1 and 3=0. 

w 


Fig. 4. 



Fig. 5. 


Figure 5 displays the rise in criti- 
cal angle x* with increasing Mach number 
and with increase in 8 from 0 to 0.2. 

Further increase in 6 has virtually no 
effect on the value of x* • 

The next series of calculations was /85 
performed for a maximally cooled surface 

at T = T = -50°C {T is the static tem- 

w e e 

perature, i.e. , 

This situation corresponds to some 
extent to a hypothetically possible ini- 
tial stage of aircraft flights in the 
atmosphere (when the aircraft has not 
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had time to heat up) . The results 
of the calculations are shown in 
Fig. 6. The sharp increase in the 
conditional transition Reynolds 
number can be seen with increasing 
Mach number (for Re t > 10 for 
M >_ 3) . Here the effect of 3 on 
Re^_ and the position of the neutral 
curve weakens. For M > 2.3 the 
frequency determining the transition 
rises with increasing 3 . As research 
showed, beginning at M = 1.5 the 
pressure gradient in the 3 = 0-0.4 
range has no effect on the value of 
the critical angle x*- 

Additional calculations at 


Fig. 6. 

T(w) = T e = -50°C ; T w = 

- 3 = 0; 

— 0.2; - • - 0.4. 


Re^iO' 2 '- 

2 


0 

0,6 0,8 1.0 



M = 4 and different 3 and T estab- 

w 

lished a moderate increase in Re t 

t 

with decrease in T to a certain 

w 

value, but with stronger cooling 
there is an abrupt increase in 
the transition Reynolds number. 

Figure 7 shows a comparison between 
the transition location calculated from 
an d "wind tunnel" data [13] with 
change in the temperature factor. 


Fig. 7. 

M = 4, 3 = 0, Tq = 300K; 
- e 2 - ' , o [13] . 


Thus, the possibility of estimating 
the position of the transition using the 
e n method (n - 9 for flight tests on a 
cone) has been shown by comparison with 


experimental data (Figs.: 4 and 7). It has been shown that while for 


3=0 and T =1 the relation between Re^ and M is determined by 
w t 

the dependence of both maximum coefficients of increase and critical 
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Reynolds number on Mach number, for 3 > 0.1 the location of the trans- 
ition is determined chiefly by the coefficients of increase. 

It has been found that whereas for a heat-insulated surface 
the pressure gradient exerts a substantial effect on the flow Stability 
characteristics and the location of the transition, with strong surface 
cooling the effect of pressure gradient diminishes. The significant 
stabilizing effect of the temperature factor on the transition location 
has been confirmed. The value of critical angle x* in the 3 = 0.2-0. 4 
range is virtually unchanged. 

The author is grateful to A. A. Maslov for his assistance in 
this work. 
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